t 


) 


4  V 


I 


1  ■ 

L  IK 


HIM 


i 


I 


r 


I  ■ 
I 


4 

1 


ilu  ■ 


'*  > 


■ 


■ 


■  * 


I  | 


H I 


I  rV" 

*  -  I  IH  1 1  I   I  H 


r 


■  * 


1 


I 


BOSTON  UNIVERSITY 
GRADUATE  SCHOOL 


Thesis 

NUMERICAL  INTEGRATION 
(Not  including  applications  to  differential  equations 
and  related  topics.) 
by 

Marion  Morse 

(B.S.  in  Ed.,  State  TeacherTs  College,  Bridgewater) 
submitted  in  partial  fulfilment  of  the 
requirements  for  the  degree  of 
Master  of  Arts 
1944 


f 


Approved 

t>7 
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y 


Second  Reader  . .  (J  JUltfA.™  •  ••^  •  *w  

Professor  of  G$/\^<yu^ 


c  .    Solution  by  Neddie's  rule. 


w  =  k 


30.24 

D  =  ydt  Z    TTT  (5. 7+15+6. 21+27+5. 3^2.32*5.30+26*6. 20+15+5. 8) 


100 


(773)  =  7.73  mi 


Example  2.  ZJve^luate 

dx 
x 


f 


3 

4 

5 

6 

7 

8 

9 

y 

.111 

.062 

.040 

.028 

.020 

.016 

.012 

a.    Solution  by  tra ^ezoidal  formula. 


f  2*  Z  i  (. nij..  1244.080+.056+. 040*. 032+. 012) 
=  |  ( .455)  -  .228 

b.  Solution  by  Simpson's  one-third  rule. 
 ■  

(  %z~  flllH(  .062+. 028*. 016)+2(.040+. 020)4.012] 

h  x      1  L 

=  -(.667)  1M 

c.  Solution  by  R  eddle's  rule. 

j    ^f-  ;-y0  r.lll+5(.002)v.042*.6(.02P)|..020f5(.016)t-.012j 

Z  _L  (.743)  z  .223 
10 


r6 

-jcample  3.    Find  the  approximate  velue  of  J     ydx  given  the  following  table 

w  here    y  -  V/o ♦  * *■ 


X 

0 

1 

2 

3 

4 

5 

6 

y 

3.162 

3.317 

3.742 

4.359 

5.099 

5.916 

6.782 

a .    Solution  by  trapezoidal  rule. 

I     ydx  ~  (1.581+3.317+3.742+4.359+5.099^.91643.391) 
-  27.405 


b.    Solution  by  Simpson's  one-third  rule. 


(  6ydx  =-    [3. 16244(3. 317+4. 359+5. 916)4-2(3.742+5.099) +6. 782) 


-  ±  (81.994)  =  27.331 
3  — — — 


• 


c .    Solution  by  Reddle's  rule. 

j  ydx  =  ~    jj3. 162+5(3. 317)+3. 742+6(4. 359)+5. 099+5(5. 916) +6. 782] 
=  v|  (91.104)  -  27.331 
Example  4.     In  travelling  by  automobile  a  passenger  checked  the  speedometer 
mileage  by  reading  the  rate  of  speed  in  miles  per  hour  every  10  minutes. 
The  readings  were 


Time 

8:30 

8:40 

8:50 

9:00 

9:10 

9:20 

9:30 

9:40 

9:50 

10:00 

10:10 

10:20 

10:3) 

Speed 

0 

10 

12 

15 

22 

24 

30 

38 

32 

25 

22 

10 

0 

Estimate  the  distance  travelled  by  car  between  9:00  and  10:00  ofclock. 

Assume  that  the  given  ordinates  belong  to  a  continuous  curve. 

a .    solution  by  trapezoidal  rale. 
-  1 

n  =  6,  w-  -g- 


/ 


10 


Sdt  =—  (15+44*48+60*75*64+25) 


=  ^2  (332)  "  27.7mi 


b.    Solution  by  Simpson's  one- third  rule. 


n  =  6,  w  =  - 


Sdt  =  —  (15+88+120+126+48+76*25) 

9  ie 


_  JL 
18 


(500)  Z  27.7  mi, 


c.    Solution  by  Reddle's  rule. 
1 

n  =  6,  w  :  — 


10  " 
Sdt  =  j±  (15+110+24+180+38+160+25) 


-  ^  (552)  =  27.6  mi. 

Example  5.    The  length  of  an  "indicator  card"  is  2.4  inches.    The  width  of 

the  diagram  at  intervals  0.2  inches  apart  are  0,0.20,  0.42,  0.53,  0.62, 

0.83,  0.89,  0.99,  1.00,  1.00,  0.98,  0.87,  0.  Determine  the  area  of  the 

"indicator  card"  and  divide  by  length  of  card  to  find  the  "mean  effective 
pressure'.'  (Dia<jwm)pa<jea3a.) 


22. 


a .  Solution  by  trapezoidal  rule. 

n  -  12,  w  =  0.2 
/-2.4 

ydx  =  .2( .20+. 42+. 53+. 62*. 83*. 89+. 99+1. OOfl. 00+. 98+. 87) 

J0 

-  .2(8.33)  =  1.67 

.7  mean  effective  pressure 

b.  Solution  by  Simpson's  one -third  ru3e  . 
n  =  12,  w  -  0.2 

f 2*tdx  _  ^  j   .80+2. 3.2+3. 32+4. 06+4. 00+3. 48+. 84+-1. 24+1. 78+2. 00+1.90 
^0 

Z  —  (25.60)  =  1.71 
3 

.7  mean  effective  ^pressure 

c.     Solution  by  Meddle  's  rule. 
n:12,w;0.2 

J ,2.4 
ydx  =  .06(1.00+0.42+3.18*0.62*4.15+1.78+4.95+1.00+6.00+0,98 
b 

+  4.35) 

Z  .06  (28.43)   ~  1.71 

.7  mean  effective  pressure 

Example  6.    The  ordinetes  of  the  boundary  of  the  deck  of  a  ship  are  0,  6,28, 

34,40,  42,  44,  42,  38,  36,  28,  6,  0  respectively,  and  the  common  distance 

between  them  is  21  feet.    Find  the  area  of  the  deck.  (Di^r*^  p<acj(:. 

a .    Solution  by  trapezoidal  rule. 

n  =  12,  w  r  21  ft. 
f252 

ydx  =  21(t6+28+34+40+42+44+42+38+36+28+6) 

'0 

=  21.344  r  7224  sq.  ft. 

.-.area  =  2(7224)  =  14,448  sq.  ft. 


23a 


b.    Solution  by  Simpson fs  one-third  rule. 

n  =  12,  w  r  21  ft. 

252 


[       ydx  =  Hi  ^(6+34*42+42+36+6)42(28^0+44*381-28)] 

Jo  SL 

-  7  1  (4)(166)V(2) (178) J  =  7140  sq.  ft. 

area  =  14,280  sq.  ft. 

Solution  by^eddle's  rule, 
n  Z  12,  w  r  21  ft. 

ydx  =  -gi  (5.6+28+6.34+40+5.42+2.44+5.42+38+6.36+28+5.6) 
=  6. 3(30+-28+204+40f210+88+210+38f 216+28+30) 

-  6.3(1122)  =  7068.6  Sq.  ft. 


area  =  14,137  .  sq.  f t. 
5! 


Example  7.  Evaluate 


J51  x 


a.    Solution  by  trapezoidal  formula. 

■55 

S  =  .009803922  * 
x  .019230769 
.018867925 
.018518519 
.009090909 


J5l 


.075512044 


b.    Solution  by  Simpson's  one-third  rule. 


dx 


I    —  =  .019607843  =  .019607843 

'51  X  4( .019230769)  =  .076923076 

2(. 018867925)  =  .037735850 

4(. 018518519)  =  .074074076 

.018181818  =  .018181818 


i  (.226522663) 
3 

=  .075507554 


The  correct  value  is  .075507552 


*Valuee  found  in  "Barlow's  Tables." 


• 


12.    central  Difference  Integration  Formula. 

a.    A  rapidly  converging  integration  formula  can  be  derived  by  inte- 
grating Stirling's  interpole tion  formula  (8:1)  from  x  r  Xq-w  to  x  =  Xgtw 

or  from  u  =  -1  to  u  =  1  where  u  =  x~*0 

w 


Then 


f  ydx  =  w  J  (y+uAy.-i+Ayp  A2  f  u(u2-i)  a3y_2i  A3y.i 
J-w        -dj  2  m  °  -l  3:  

u2(u2-l)  A4            u(u2-l)(u2-4)  /  5y  tA5y 
— 4?       a  y_2  -t-  51  ^  -2_  =£_ 


u2(u2-l)(u2-4) 

^6y_3+ 

3 


61 


2  3  4      2  3J2 


.3    <K  „6      ^,4  «  ,5. 


M»    .  u  )AV^    Mu_  _  5u*  +  4u^'y_3-^y_2 
4! 


5  J 12 


n7       5u5      4U3     £  y  T  1 

:  w^yc^y^  ^  (f  §  )  A4y.2t  &  |  -it  |)  A6y_3f  -J 

=  2w  fy0  f  7&2y      -    -i-£4y_2  v  — ~ A%      t  *•*  1 
Uu      6        -1        180  2      1512"  y-3  J 

To  integrate  from  x  r    xQ  to  x  r  Xq+2w  advance  y0  to  y^,  y_j_  to  yQ, 
y-2  t0  y-l 

This  procedure  can  be  continued  until 

>C2to  z  2w  ^n-x  *  *  A2vz "  ^ A4yn-3  *■  A6ynJ 

Then  by  adding  these  integrals  we  obtain 

J-  xQ*-nw 
ydx  =  2w£yi  v  y3  v  y5  f      ^  yn-i  +  |(A2yofA2y2*-A2y4* * *A2yn-2 

"i8o  ^y-i^4yi^4y3+*--  +  A4yn.3) 

(12:1)  I 

*  USS  (A6y.2^6y0+AV*--^6yn-4lI 
where  n  is  an  e-ven  number. 


b.    The  approximate  area  may  also  be  fo\ind  by  integrating  Bessel*  s 
interpolation  formula  (8:5)  from  x  -  xQ  to  x  s  Xq  +  w  or  from  v  =  ~£  to 

▼  8  i 

Thus  /X°+W  ydx  =  w/*  ISL±1±    ♦  v4  y0  +  (**  "  j>  A**»l  -  ^~ 


-  w 


v(tj  -  i)  43y.i  tvs  -  i)(vE  -  I)  a4y.2  -  a4y_, 
*  Ti  +  4i  g  * 

„  ^  -  i)  (**  -  j)  A-g  ,  (*2  -  i)  (v2  -  |)  (v2  -  |&>  A-3  -  A-b 

5!                                                 61  2 
f  v(^-i)(T2-|)(v2-|^)   ^  )dT 

3  J 


V5        6  16  41  2  6        8    T  32  5J 


7_£      259V3      225y.  46y-3  »    A67.8    +  ^  259v4  225^ 

7    ~   4      f    48  64  '  61  2  8    ~  24  64      ~  128 


4^3 

7j 


:  w 


— 2  r2  S  *  72"0"  5 


_  191  A6y-3  t  A6y-2 
60480  2 


Then 


xQ  f  2w 


ydx  s  t? 


x0  f  w 


yx  +  y2     1  A*-y0  ♦  A27i      11  ^V-i  +  A4y0 

 S-*    -  TS   o   +  T20   5  


12 

191  46y-2  ♦  46y-i 


60480 


ktai  l"1        12  2        720  2 


i9i  4 6y-i  f  <a  yo 

60480  2 


This  procedure  may  be  continued  until 


xQ  +■  nw 


ydx  *  w 
f  n-lw 


15" 


191    4  6yn-4  t  *  6yn-3 


720" 


60480  2 
When  added  these  integrals  give 


(12:2) 


x    f  mw 
o  _ 

ydx    z  w 


£*7xt,8t  W^  -^(^ 


*  aSt    *  *    A4v       *  d^n-g)  _      191  f^Cy-3  t  46y 

*  A^y0  ♦  *  •  •  ♦   A  yn-3  *  —7  '      604080    a  y- 


6 


•   •  • 


*  aS 


n-4  +   "  ; 


+  A   7-1  ^ 

n 

where  n  can  be  even  or  odd. 

This  can  be  simplified  by  changing  the  even  differences  to  odd 
differences  in  the  following  way: 


A  y-1  =  A  y  -  A  y_: 

a  2yo  *  ^  yi  -  4  yc 


4  yn-l=  47n-  4  yn.! 

zi4y.2  =  ^3y-i  -  A  3y-2 


d4yn-2  =  43yn-i  -  A3yn-2 
A6y.3  =  d5y.2  -  A5y-3 


a6yn-3  =  /:5yn-2  -  ^5yn-3 


and  substituting  in  the  preceding  formula  (12:2)  as  follows: 


'xQ  +  mw 
ydx  -  w 


y0 


11  (  ^3y.2  *  43y-i?  +  i9i  (^5y.5  t  /y-2  m 


•)  * 


720               2  60430  2  12 

11    ,  43yn-2  ♦  ifl^n-l  5V    „  *  A  5V  c 

+  W~<  i  )  -  JSkr^  7n-3  g  *  yn"g) 


60480 


xQ  *  nw 


r 


or  (12:3)  /      ydx     =  w 


f  *  yx  *  y2  *  •  .  .  ♦  yn-i  +  ^4^y»-y4y» 
^'y.i  +  ^yP      ii  /yn-2  +  43yn-i  _  /y-2  *  aV-i 

2         '       720  2  2 

191  ^5yn-3  *•  45yn-2    45y-3  +•  /y-2x 


"  60480 ;  2  2 

In  (12:1)  the  graph  of  the  function  is  replaced  by  ^  arcs  of  parabolas 

of  sixth  degree.    In  (12:3)  the  graph  of  the  function  is  replaced  by  n  arcs 

of  sixth  degree  parabolas. 

If  in  formula  (12:1)  the  second  differences  are  replaced  by  values  in 
terms  of  y  it  will  be  seen  that  the  terms  containing  fourth  and  sixth 
differences  are  correction  terms  for  Simpson's  one-third  rule  (11:3). and 

For:  2wjryi  *  y3  *  .  .  ,  +yn-1  +  ^U^q  *  42y2  d2yn-2) 
would  become  2w  [yx  f  y5  ♦  ♦  •  .  4  yn.1  *  -i(y2  -  2yx  *  yQ  *  y4  -  2y3  +  y2  + 

.  .  .  *  yn-2  -  2yn.!  ♦  ynj 
y0  *  4(yi  t  y3  *  .  .  .  +  yn-1)  t  2(y2  +  y4  *  .  .  .  *yn.2)  t>  yn 
which  is  (11:3) 

Formula  (12:3)  represents  the  trapezoidal  rule  with  correction  terms. 


=  3 


c 


12.  Examples 


X 

y 

A  y 

/  y 

A  i-7  v 
03.7 

48 

.02083333 

49 

. 02040816 

-42517 

50 

.02000000 

-40816 

1701 

51 

.01960784 

-39216 

1600 

-101 

52 

.01923077 

-37707 

1509 

-  91 

10 

53 

.01886792 

-36285 

1422 

3  67 

4 

-6 

54 

.01851852 

-34940 

1345 

-  77 

10 

6 

12 

55 

.01818182 

-33670 

1270 

-  75 

2 

-8 

-14 

56 

.01785714 

-32468 

1202 

-  .68 

7 

5 

13 

57 

.01754386 

-31328 

1140 

■  63 

6 

-1 

-  6 

58 

.01724138 

-30248 

1080 

-  60 

2  • 

-4 

-  3 

a.    Solution  by  Stirling's  formula  (12:1) 

r  55 

/       il  s  2 (.03774929}  =     +  .07549858 


t  ^.(.00002692) 


+  .00000897 


f  .07550755 


b.    Solution  by  Bessel's  formula  (12:3) 
55 

—  =       .07551204    =  .07551204 
51  X 

-  jfit 00005392)    :      -  .00000449 


.07550755 


Example  2.    For  data  refer  (11    ,  example  4) 


77 


— 4~ 

a  y 


"5" 


~~ 6~~ 

<o  y 


a  y 


A2y 


7 

4  y 


8:30 
8:40 
8:50 
9:00 
9:10 
9:20 
9:30 
9:40 
9:50 
10:00 
10:10 
10:20 
10:30 


0 
10 
12 
15 
22 
24 
30 
38 
J32 
25 
22 
10 

0 


10 
2 
3 
7 
2 
6 
8 
■  6 
»  7 
-  3 
■12 
•10 


-  8 

1 
4 

-  5 
4 
2 

-14 

-  1 

+  4 

-  9 
*  2 


9 

_  3 
-9 
9 
-2 
-16 
+  13 
*  5 
-15 
til 


-  6 
-12 

18 
-11 
-14 
+29 

-  8 
-18 
+24 


-  6 
+30 
-29 

-  3 
+43 
-37 
-10 
+42 


+36 
-59 
+  26 
+  46 


-SO 
+27 
+52 


-  95 

h  85 
+  20 
-126 
+  107 
*  25 


+180 

-  65 
-146 
+233 

-  82 


a.    Solution  by  Stirling's  formula 


10 
Sdt 


n  =  6,  w  -  -i  nr. 


.  I  [22  +  30  *  32  +  i(-5  +  2  -  1)  -  ^(18  -  14  - 
"  3  I,  6  160 

+  _L-_(-59  +  46  +  27)  =  27.78  miles 
1512   

b.     Solution  by  Eessel's  formula 


8) 


n  =  6      n  z  —  hr. 


10 


S  dt  =  i    166  -  _i(-10)  +  _11(-1)  -     151  (a) 

6  L  12  720  60480 

-  1(266  +  .833  -  .015  -  .012)  =  27. 8C  miles 
6   — — — 


31. 

• 

1 

13.    Zuler-I/aclaurin  Formula.    Let  f  (x)  r         where  v  is  independent  of  x 
and  substitute  in  the  trapezoidal  formula  (11:2).  Then 

(13:1)  if0  ev^dx  a  J  e^o  f  ev(x0  i  w)  f  ev(x0  ^)  f 

^0 

qv(x0  f  n-lw)  +  t  Qv(x0  ♦  nw)  +  ?  where  T  ^  ^  correction 

term  in  the  trapezoidal  formula. 

/x0  t  nw         v(x    f  nw)        vx          vx0  vnw 
irtr.tl'M  1/     «vx  riv  -  e      o               e    °  nr  e      (e       -  i) 
w/_                            vw  vw 

Thuo  «VX°(eVn*  "        -  t  «VX°  +  ev<xo  +  w)  4.  «v^c  *  2w)  +  .  .  .   .  + 
vw 

ev(x0  4  n^lw)  ^  ^v(x0  +  nw)  ^  t 
s  e™o  f  1  +  evw  +  e2vw  +  _  ^                     (vnw  _  ^  ? 

_  e™o(e™w  -  1)  .  ..vxo^vnw  *  tj  4  t 
ev*  -  1 

Multiplying  hoth  sides  by   

evx0(evnw  _  1} 

(13:2)         1  5  -3H         +  ™  *  VVI 

1                        evw  .  1       2      evx0(evnw  „  x) 

vw         ,      vw  .  BXv2w2      B3vV     B5yV  B7v8w8 
If  expanded  in  series               =  1  -  -  -  f    gy       -      4j      4    ^  ^ 

e     -  1 

(13:3) 

4  .  .  •  •  • 

hen  (13:3)  is  substituted  in  (13:?) 

vw      B!V2w2      B3v4w4      B5v6we  B7vew8 

1  =  1  "    2  f      2!            4!      i      61      "      e:  f 

(  vw  (  vwT 

n         vx„ ,  vnw      -  , 
2       e    o(e       -  1) 

1 

Dividing  through  by  vw  and  solving  for  T  the  result  is 

„                        Bi  vw      B^vV      B5v5w5  B7v7w7 
T  =  e^ojo™  -  i)(-    1      +    3         -  ~2         f  ) 

?,!            4:              6!  8i 

Then  f(x)  =         then  f » (x)  =  vevx,  f"(x)  =  vSevx,  .... 

so  that  T  =  22-    f(x0  +  nw)  -  ffxW  -i  JU  f  '  "  (xQ  +  nw)  -  f  "(*(,)/ 
2!     L  J         4  •    I—  J 

.5  r  ,      B  «7r- 


B5W 

e:  l 


fv(x0  t  nw)  -  fv(x0)|  4  £2LffVii(x0«w)-fvli(x0)J.  .  . 


p»se 


8!  L 

Thus  when  f(x)  =  e      formula  (11: 2;^  plus  T  becomes 

I  f(x0)  f  f(xQ  +  w)  *  f (  x0  +  2w)  +  .  .  .  t  f(xQ  *  n^lw) 
B3W 

a 


fx    +  nw 

f \x)  dx  =  w 


*  if(x0  +  nw)  -  —  ^f»(x0  t  nw)  -  f •  (x0)j  *  -frr*  (xo*n*) 


-f"  (x0) 


fv(xD  *  nw)  -  fv(xfi)    t  .  .  .  . 


ij  ♦  .  k  i  k  J 


=  w  Uy0  t  yx  t  y2  *  .  ,  .  4-  yn.1  4  tyn  -  —  (y£  *  y0) 

»  fff^  -  ,"0)  -  !g!  (,J  -  #  f  ^tfM")  •  •  •  •] 

The  B's  are  known  as  Bernoullian  numbers  and  have  the  values  B^  =  -jr  » 


B<z  ■  — i,  B r-  s    ^* »  Bn  *  — ^  «  •  .  •  • 

0      30*     5      42      7  30 


x0  f  nw 


(13:4)  Thus/     f(x)  dx  =  w 


yo 

i 


♦  >vi +  r  '  |r(y» "  y°} 


-3  »5      ,  v       v>         w7      /  vii  vii 

t  £o(y£'  -  y"0)  -  ^T0(yn  -  y0)  *  i^960^(y  »  "  y  o)#  ' 


This  is  known  as  the  Euler-lilaclaurin  formula. 

It  is  necessary  to  show  that  this  formula  applies  not  only  when 
f  (x)  -  e^  but  when  f  (x)  is  any  arbitrary  function.    To  do  this  assume  that 

V«iX 

it  is  possible  for  the  sum  A]_ev^x  +  Age  2  where  A^f  Ag,  •  .  • 

and  v^.,  v      .  ,  .  .  are  constants  properly  chosen  to  represent  a  close 
approximation  to  f (x)  between  finite  bounds  x  «  rQ  and  x  s  xQ  f  nw. 
A  and  v  are  independent  of  x.    Then  if  formula  (13:4)  is  applied  to  each 


4  « 


t 


t 


ViX  VpX 

term  A^e        i  Age        >  •  •  •  •  and  the  results  added  it  can  be  seen  that 
formula  (13:4)  may  be  applied  to  any  arbitrary  functionf(x) . 
Note:    The  terms  in  the  Euler-LIaclaurin  formula  beginning  with  'jsji^n  "  y'o) 
form  an  assymptotic  series,  which  is  an  infinite  series  which  converges  for 
a  certain  number  of  terms  and  then  begins  to  diverge.    Thus  in  computing 
with  this  formula  it  is  important  to  know  what  term  to  stop  with  in  order 
to  get  the  most  accurate  result.    One  should  stop  with  the  term  just  before 
the  smallest. 
13,  Zxair.pl e s 

f  55 

U    Svaluate  /  — 

a.    Solution  by  Euler-IJaclaurin  formula 
y  :  I        wsl  n  =  4 

X 


■        2  51      52      53      54      2  SB1      u<     652      512>  +  ^  C-  6g4  *  ^71 


•00980392 
.01923077 
.01886792 
.01851852 
.00909091 

.07551204  =  +  .07551204 
-  .00033058 
+  .00038447 

-  _1_(+  .00005389)-  -  .00000449 

12   

+  .07550755 

loge  55  =  4.007333185 
■logQ  51  5  3.931825633 


0.075507552    correct  value 


*  The  reason  for  this  is  given  by  Cherlier's  rule  in  Chapter  IV.  Article  25. 


34. 

f*,  r 

Example  2.                Evaluate    /    ylO  f  xc  dx 

a«    Solution  by  Euler-Maclaurin 

y  r  VlO  4  X2                    W  -  1        n  =  6           f  *{x)  m          *  , 

f 6                                           r-  fx>  -  3x3 

/    y  dx  =    1.581139                                            1           *  /-  ZTT 

JO                3.316625  7(10+x2)5 
3.741657 
4.358899 
5.099020 
5.916080 
3,391165 
27.404585    =  27.404585 
-  _i  (     .841173)  »  -  .070097 
12  27.33448e 

Example  3.    Evaluate  J^f  =  ^  ik  +  ijo  +  W  f  2^0  +  2*0  *  25^ 

W=10                                  ,  10(      1      ♦      1  ) 

r  12             p  2 
n  :    5                                           "(250}^  (200) 

=    .2232185  *          000009)  r  .2232260 

14.    Gregory* s  Formula.    The  Suler-Maclaurin  formula  may  be  expressed  in 

terms  of  differences  instead  of  derivatives.    Formulas  for  the  derivatives 

of  y    may  be  found  by  differentiating  Newton* s  interpolation  formula  (7:1) 

and  putting  x  a  0  in  the  set  of  equations  obtained.    Those  for  f(n)  may 

be  obtained  by  forming  a  difference  table  and  reversing  the  order  of 

1 

differences  so  that  the  last  is  first,  next  to  the  last  second,  etc.,  and 
then  applying  equations  for  derivatives  of  f(0). 

Some  of  the  equations  thus  obtained  are  as  follows: 

3         ,    4         i  .5 


^0  = 

4y0  -  \  A\  » 

= 

*yn-l  +  I4'7* 

-  |^4y0  ■ 

3  „, 
w  y-r 

A3yn_3  *  1  ^y 

5  v 
Wyn  = 

A5y0 

5  v 

A  y3  = 

5 

A  7n-5  +  '  •  • 

etc. 

•    •    •  • 


1  t\  

.  7  ,5„ 


Then  these  are  substituted  in  the  Euler-l&claurin  formula  it  becomes 

Vo  .  y, 


x0  *  nw 


(14:2)    /       y  dx    =  w 


x. 


-       *  V*  *  A>  -  ^(  4 5yn.5  -  ■ 
fi^Ve  *  A0)  ] 


24192 

This  is  known  as  Gregory's  formula  and  has  the  advantage  of  not 
requiring  ordinates  beyond  the  limits  of  yQ  and  yn  in  order  to  form  differ- 
ences. 


f 


36. 


14.  Examples. 


55 


/  — 
51 


A7 

3 

z 

y 

A 

4  y 

51 

.01960784 

-37707 

52 

.01922077 

-36285 

1422 

-77 

53 

.01886792 

-34940 

1345 

-75 

54 

.01851852 

-33670 

1270 

55 

,01818182 

102  T  52      53      54  110 


12 


(.00004037)  - 


00002  692) 

24 


 i£(.  00000002)7 

720  _J 


.07551204 
-  .00000336 

.07550868 


Example  2.       Use  data  from  (11  example  6) 
Solution  by  Gregory's  formula.       (See  next  page.) 


z 

y 

dy 

A 

A  y 

*  6 

A  y 

0 

0 

6 

1 

6 

22 

16 

-32 

2 

28 

6 

-16 

+16 

+48 

-68 

3 

34 

6 

0 

-  4 

-20 

+28 

+96 

4 

40 

2 

-  4 

f  4 

+  8 

-16 

-44 

5 

42 

2 

0 

-  4 

-  8 

+  14 

+30 

6 

44 

-  2 

-  4 

+  2 

+  6 

-  4 

-18 

7 

42 

-  4 

-  2 

+  4 

+  2 

-14 

-10 

S 

38 

-  2 

+  2 

-  8 

-12 

+32 

+46 

9 

36 

-  8 

-  6 

+12 

+20 

-42 

-74 

10 

28 

-  2 

+  6 

-10 

-22 

11 

6 

-  6 

-  4 

12 

0 

=  12 

f  252 

y  a* 

t/0 


w  .  21 


2l[344  -  M-  6  -  6)  -  4  +  16)  -  -  10  +  32) 

_  -JL(-  22  +  48)  -  — -  42  +  68)  -      275 f     74  +  96). 
160  60480 V  '      24192 1  f 

21(344  +1      1  39      11219  3025) 

1  2  "  360  "  80  "  30240  "  12096 

21  [344  +  1  -  .50  -  .58  -  .49  -  .37  -  .25^ 
21  (342.81)  ■     7199.01  sr^.  ft. 
Area»14398.02  sq.  ft. 


Example  3.    One  side  of  a  field  is  straight  and  of  length  200  feet.  At 
distances  increasing  by  20  feet  from  one  end,  the  width  is  measured  and 
found  to  he  80,  110,  110,  140,  147,  160,  155,  110,  110,  80,  50  feet. 
Find  the  approximate  area  of  the  field. 


X 

y 

ay 

V 

43y 

44y 

A  5y 

46y 

0 

80 

20 

110 

30 

40 

110 

0 

-30 

60 

140 

30 

30 

60 

80 

147 

7 

-23 

-53 

-113 

100 

160 

13 

6 

+29 

+  82 

+195 

120 

155 

-  5 

-18 

-24 

-  53 

-135 

-330 

140 

110 

-45 

-40 

-22 

+  2 

+  55 

+190 

160 

110 

0 

+45 

+e5 

-107 

-109 

-164 

180 

80 

-30 

-30 

-75 

-160 

-  53 

+  56 

200 

50 

-30 

0 

+30 

+105 

+265 

+318 

w  z  20,  n  =  10 

^200  r 

JQ  y  dx  =  20  [  1187  -  60)  -  -A<-  30)  -  30>  ~  ^  *)  " 

.  __§£3(70)  -  _275(18)7 
60480  24192  J 

=  20  [ll87  f  5  +  1.25  +  .79  +  .15  -  1.00  -  .20 j 
=  23859.80  sq.  ft. 

15.    Woolhouse's  Formula.    The  formula  is  used  to  find  the  value  of  the 

sum  of  a  large  number  of  terms.     If  the  interval  chosen  between  the 

successive  values  of  the  argument  in  the  Euler-Maclaurin  formula  is  w  and 

it  is  desired  to  reduce  this  interval  to  —  then  the  formula  (13:4)  becomes 

m 


(15:1)  £ 


'x0  +  nw 


When  formula  (15:1)  is  subtracted  from  m  times  the  formula  (13:4) 


the  following  result  is  obtained 

(15:2)     y    +  y       +  y       +  .  .  .  f  jr        +  y    z  m(y    +  jx  f  ....  4y  ) 

1/m       2/m  >ili 

m 

2  4  3 

-  m  "  1{j    +■  y  )  -  (m    -  1)W(T«  -  y»)  x  (m    -  l)w  ,  w  _        )  + 
2     vyo      V  12m         n       o'  r      720  mA    iyn     yo  ,T  " 

Thus  by  increasing  the  number  of  intervals  so  that  the  new  intervals 

are  —  times  as  great  the  sum  of  a  large  number  of  values  can  be  found  from 
m 

a  small  number  of  values. 

16.    Lubbock's  Formula.     If  differences  (14:1)  are  used  in  place  of  derived 
functions  in  (15:2)  the  following  formula,  called  Lubbock's  formula,  is 
obtained 

(16:1)    y0  +  y^  f  yg/m  +  .  .  .  +  y^f  yn  =  m(y0  +  yx  +  .  .  .  +  yn) 


39. 


m 


720  m 


-  ^-1"9m'-1)  (  * \-4  ♦  A„)  -  •  • 


480  me 


Examples. 


250 


1.    Determine  value  of 

n  =  200  u 

a.    Solution  by  Woolhouse's  formula 
m  s  10  w  =  10 

250 


I  -  10 (-1-  T  -L-  f  -i-  +  -JL  t         +  -i-)  -  9(  1    +  JL) 

n  ^200  n  200      210      220      230      240      250 ;      2*200  +  250 ' 


99  (  _  1 


t^~) 

2502  goo2 


♦The  numerical  coefficients  ^^iffi*  listed  as  An(r)  *nere  r=m  are  found 

in  "Tables  of  the  Higher  Mathematical  Functions"  by  Harold  Davis, Table  34. 


s  .005000000 

.004761905 

.004545455 

.004347826 

.004166667 

.004000000 
10 (.026821853)  =  +  .26821853 

.005000000 


.004000000 

-  i  (.009000000; 

-  .000016000 


s  -  .04050000 


99 


f  .000025000 


"  ff  (*  .000009000)  =  -  .00007425 

t  .22764428 

.22764430  was  the  correct  value  found  when  the  reciprocals  were  added. 
1.  b.  Solution  by  Lubbock* s  formula 
m  z  10         w  :  10 

250  Q  QQ 

£  1  z  io(y0  t  yx  +  •  •  .  *  y5)  -  |(  y0  +  y5)  -  1|§(^y4  -4y0) 

n  r  200 


-_H(^y3,^0)-i|lie|9|W3yg.A0, 


y 

ay 

A27 

A  y 

a  y 

.005000000 

-  ,000238095 

.004761905 

-  .000216450 

.000021645 

-  .000002824 

.004545455 

-  .000197629 

.000018821 

-  .000002351 

+  .000000473 

.004347886 

-  .000181159 

.000016470 

-  .000001978 

+  .000000373 

.004166667 

-  .000166667 

+ 

.000014492 

.004000000 

250 


A  =  0.26821853  -  |(  .009000000)  -  -—(+. 000071428)  -  _?!(+. 000036137 


n  =  200 


240 


-  Ml899l(  +  .000000846)  -  ||^i( . 000000846 )  =  .22764234 


480000 


80 

2.    Determine  value  of  (1  +  x2) 

x  s  50 

a.    Solution  by  Woolnouse's  formula 


x  z 


80 
50 


12 


I  - 


80 

To" 


m  s  10      w  -  10 

=  10(2501  +  3601  f  4901  +  6401)  -  |(2501  *  6401) 
-  174040  -  40059  -  495  =  133486 

b.    Solution  by  Lubbock's  formula. 
=  10(yo  i  yx  +  y2  +  y3)  -  |(  y0  +  y3)  -  ^(dy,,  -dy0) 


(100  -  160) 


-  2iu2yif  Ao} 


y 

a  y 

A2y 

2501 

t  lioo 

3601 

+  200 

t  1300 

4901 

♦  200 

t  1500 

6401 

£    z  174040  -  40059  -  330  -  165    z  +  133486 


80 
i  =  50 
Example  3. 

Evaluate 


250 


—   by  Euler  Maclaurin  formula  and  Woolhouse  formula 
x 


J  200 

m  z  10,  w  s  10 

yQ  ^  y  .    +  .  .  .  +  y    =  .22764428  by  (16  example  1) 
i/m  n 


Then      H  Ao  +  »     /*5°ta  .1,1    .    1  ,       1,  _1     +       1  , 


•  .22764428  -    -|  (+.009000)  -  -~(+. 000009) 


=  ,22764428  -  .00450000  -  .00000075  a  .22314353 


41. 


It  is  noted  from  (13  example  3)  that  when  n  s  5  and  w  s  10  the 
answer  was  .2232860 

17. TKe  Newton-Cotes  Formula.    If  it  is  desired  to  find  the  integral  between 

the  limits  x  s  a  and  x  a  b  the  limits  may  be  changed  to  -1  and  +1  by  the 

(+1 

transformation  x  =  a  *      -  ^  "  a  £  •    To  integrate  the  function/  f(x)  dx 


22^  y-i 

try  to  find  a  sum  of  terms,  say  HQf(hQ)  +  H.f(h.)  +  H_f(h  )  with 

11  Ai        ^  n 

H0,  H^,  Hgl  •  •  •  ,  Hn  constants,  which  will  represent J  f (x)  dx  as  nearly 
as  possible.    If  the  sum  is  to  be  equal  to  the  integral  exactly  as  long  as 


f (x)  is  any  polynomial  of  degree  n  at  most  then  f(x)  may  be  represented  by 
Lagrange's  formula  (9) 

3      Ffx)  f<hk> 
»     *>      ,V  \  *  „tt  \  '    711611  ^  integrating 
k  s"o  1        h*'  F  <hk' 

ffl  3      f<hk>    f*1    F(x)  3 

.'.  Ev  .  -4         /    1  -Ii*L  dx 

k     FlSJ-y^    x  .  hk 

If  the  h's  are  taken  at  equal  distances  apart  so  that  h_  =  -1,  fcu  =  -1  +  ~, 

h2  ,  -  1  f  i,  .  .  .  h    -  1 

n 

then       -HiL  =  (x  f  1)  (x  4  1  -  £)  .  .  .  fx  f  1  -  2(k  ~ 

x  -  hk  H  L  n  J 

[x  + 1  -  my  uj ...  (x  - 1) 

If  z  =  |(x  ♦  1)  then 


2 


=  (£)"f z(z  -  l)(z  -  2)  .  .  .   (z  -  k  4  l)(z  -  k  -  1)  .  .  .  (z  -  «)7 
-  n    L  J 


Now  —i- —  -  1  

(hk)      ^  -  h0)<hk  '  nl>  '  '  <nk  "  hk-l><hk  *  W  '  -(nk  "  V 


=  (^i)n"k(I)%J(a-k)l  (-l)n"k2nk!(n-k): 
Thus  J      fix)  dx  =  H0f (-1)  f  Hxf (-1  f  1)  t  .  •  •  f  V (1) 

where  (17:1)  Hk  -  J^z'1]  ■  '  *  .  .  .  (z-n)dz 

or  if  h0  s  ax,  hk  =  a  f  kw  then 
fxQ  +  nw 

(17:2)    /     f (x)  dx  .  H()f(x0)  f  Hxf (xQ  *  w)  +  Egf (xq  +  2w)  *  .  .  .  Hnf(xo+nw) 
%/  xo 

n-k     /*  n 

where  Hi.  -   A  /    z(z-l).  •  •  (z-k+1) (z-k-1)  .  .  (z-n)  dz 

*  k  -  k! (n-k) : J o 

Some  of  the  values  of  are: 
n  ,  1     EQ  =  H1  s  | 

n  =  2      H0  =  H2  =  §  5  Hx  =  £ 


n  =  3     HQ  -  H3  =  —  ;  Hx  r  H2  =  — 


14w 

n  =  4  HQ  =  H4  s  -jg 
ns5  H0,H5S|g 


n  s  6      Hrt  s  HA  =  41w 


-  w    -  £4*  •  h    -  8w 

Hl  s  H3  =  "45  »  H2  =  15 


tt  rr  125W     .     TT  TT  125W 

Hl  =  H4  ■  —ge  '  H2  =  H3  = 

54w  .  TT       „        27w  .  „  204w 


96  '     2  -    3  *  144 

°  =  "6  =  Tio  '  Hl  3  H5  ~  ^  5  Hs  =  H4  5  I5&  5  ^  =  "T05 
WKen  the  Newton-Cotes  formula  is  used  it  is  simpler  to  divide  the  interval 


from  a  to  b  into  a  number  of  sub intervals.    The  formula  is  then  applied  to 
each  subinterval  and  the  value  of  the  whole  integral  found  by  taking  the 
sum  of  the  separate  terms  obtained. 

The  trapezoidal  rule,  Simpson's  rule,  and  Wtddle's  rule  can  also  be 
obtained  from^Newton-Cotes  formula  by  taking  n-1,  nr2,nr6  respect- 
ively. 


17.  Examples. 

1.    Use  data  in  (11  example  4)  and  solve  by  Newton-Cotes'  formula, 
n  :  6,  .     i,  H0  S  Hg  .  _«         =  Hg  ;  _|,  Hg  ;  H4  :  %  =  ^ 


10 

S  dt  : 


41 

840 


15.  _£.22  ,   £.24  .  _3±.30  ,   £.38  .  __£.32  .  _ii.25 

35 


280 


105 


280 


35 


840 


s  .732  +  5.657  f  .771  f  9.714  +  1.221  +  8.228  f  1.220 
Z  27.543  miles 

"l.5  3 

2.    Evaluate  /         •   dx  by  NewtOE-Cotes  formula. 

0       ex  -  1 

This  integral  is  important  in„Debye  theory  of  heat  capacity  of  solids. 


X 

f (x)  =  x3/eX  -  1 

u 

0 

0.25 

0.055013 

0.50 

0.192687 

0.75 

0.377686 

1.00 

0.581977 

1.25 

0.764280 

1.50 

0.969357 

3x.25 


9x.25 


n  =  3,  w  =  0.25,  H0  s  H3  s  ^  «  Hg  -  - 


0.75  3 
x 


e*  -  1 


=  *  ^(0.055013)  f  ^(0.192687)  +  377686) 


r. 015472  +  .054193  +  .035408  =  .105073 
—  dx  :  -^(. 337686)  f  li^. 581977)  +  ^f5^. 784280)  f  -^(.969357 
=  .031659  +  .163556  f  .220579  +  .090877  =  »506671 


1.50  3 
x 


e*  -  1 


dx  =  .105073  f  .506671  =  .611754 


3,  Evaluate 

+  1       dx  tt  :  E,  I  :  1,  H    :  L  :  I  ,  H,  :  * 

TTZ  o       2      3        1  3 

\  z  -1,  h1  =  0,  h2  =  1 


/  1      x  +  2      3  T  3  f  9  "  9 

The  correct  answer  is  1.09861.    The  above  answer  is  ineccurate 
because  of  the  small  number  of  subdivisions. 

18.    Gauss'  Integration  Formula.      In  Gauss's  forrula  the  h's  are  chosen 
at  such  distances  apart  that  the  formula  of  Newt on- Cotes  will  have  the 
greatest  possible  accuracy.    As  there  are  2n  +  2  constants,  the  h's  and  H's 
can  be  chosen  so  that  the  formula  is  exact  when  f (x)  is  a  polynomial  of 
degree  2n  +  1  at  most.    If  F(x)  is  a  polynomial  of  n  f  1  degree  denoted  by 
(x  -  hQ)(x  -  h^) (x  -  h2)  •  •  •  (x  -  hn)  as  in  Newton -Cotes' formula  and 
<^(x)  is  a  polynomial  of  degree  c  at  most  then  it  is  necessary  that 
ftl  ftl 

/    f  (x)  dx  -  /      <p  (x)  F(x)  dx  s  0 

/VI 

As  J    (p  (x)  Pn+^(x)  dx  =  0  when  Pq+^(x)  is  a  Legendre  polynomial  and  <p{r) 

is  any  polynomial  of  degree  n  at  most  then  the  conditions  are  satisfied  if 
F(x)  -  Pn+l'x)  and  the  nuD3:ber8  n0»         •  •  •  nn  in  tne  formula  (18:1) 


(l8:l) 


}  ^  f (x)  dx  =  H0f (h0)  +  Hxf (hx)  .  .  .  Hnf (hn) 


*  The  proof  is  given  in  "A  Course  of  Modern  Analysis'*  by  S.  T.  ""hittaker 
and  G.  N,  Watson,  pg.  305  -  15:14 


* 

are  the  roots  of  Legendre  polynomials  of  degree  n  f  1.      The  H*s  are  then 

)    .  l       J(x)    dx    as  in  Newton- Cotes'  formula, 
iouna  oy  nk  _  F' (hk) J_±  x  -  h^ 

Some  of  the  values  of  h's  and  H's  are: 

Ho  =  Hl  =  1 


n  =  2         ho  =  -7T  0  S.    1  '  ? 


*o 

1 

=  '  7T 

*1 

1 

=  rr 

n 

ho 

hl 

-  o 

h2 

=  7^ 

of 

R_  and 
m 

found  from  the  relationships  Hv  -         and  h.   -  2V^-a.  k-m.  R    and  v    are  the 

A 

-alues  in  the  formula   /    f (x)  dx  =  RQvc  *  Riv1  +  *  *  ,+Rnvn  where  tGe  i^ter- 

val  is  from  0  to  1  instead  of  from  -1  to  +1. 

Causs's  formula  may  be  used  for  another  purpose.     It  may  be  used  for 

1  The  Legendre  polynomial  is  defined  by  P_(x)  =  — i—  - — (x2  -  l)n.    Some  of 

n  2nnJ  dx11 

the  values  are  F0(x)=l,  P1(x)=  x,  Pg(x)  =  £(3x2-l) ,  P3(x)=i (5x3-3x) .  The 

values  of  F0(x),  P1(x)  ,  ...  are  found  in  "Tables  of  Functions'*  by  Dr.  Eugene 
Jahnke  and  Fritz  Smde.  pg.  ICS  and  the  roots  are  found  by  setting  the  -poly- 
nomials    equal  to  zero. 

**Appendix  3 

***"The  I.^thematics  of  Physics  and  Chemistry1*  by  Henry  I/Iargeneau  and  George 
..oseley  Murphy,  pg.  463. 


finding  the  test  times  for  taking  measurements,  as  for  example,  to  find 
the  best  times  of  day  to  take  the  temperature  in  order  to  obtain  the  best 
value  for  the  mean  daily  temperature  if  a  certain  number  of  measurements 
are  to  be  taken. 

The  mean  value  of  a  function  is  given  by  the  formula 

A> 


y  dx 


'a 


ym  =  b  -  a  rb 
The  accuracy  depends  on  the  accuracy  of  /     y  dx.    The  most  accurate 

valua  is  found  by  taking  ordinates  given  by  Gauss'  formula. 

If  the  interval  x  z  a  to  x  z  b  is  transformed  to  i  :  0  to  x  :  +1  then 
X  s  a  f  (b  -  a)v. 

When  it  is  desired  to  find  the  best  times  of  day  to  take  the  measure- 
ments of  temperature  in  order  to  obtain  the  best  value  for  mean  daily 

temperature  the  formula  may  be  written  as 

T24 

/    t  (dt) 

Tw  =  J 9.   Then  t  -  24v 

m  24 

* 

The  v's  are  found  in  the  table    for  n  measurements  and  the  t's  can  then 

be  determined. 

IS.  Examples. 

1.    Evaluate/  — —r    by  Gauss'  formula. 
J-l       +  1 

n  =  2,  vQ  z  .11270166,  y%  z  .5,  vg  =  .88729833,  HQ  :  Rg  =  ]|»  Rl  :  {  J 
hk  «  2vm  -  1,        Hk  =  2  hm 

V(>o>  =  1^54  *  f        •  I  =  -5555 

*  Appendix  B 


H2f (h2)  r  2.7746  •  e 

/  41    d*    :  .4534 
/  -i  x      2  .4444 
.2002 
1.0980 

The  correct  value  is  1.09861 
Example  2.    It  is  desired  to  find  the  mean  daily  temperature  by  means  of 
5  temperature  readings.    At  what  times  ought  these  to  be  taken?    Take  the 
day  as  extending  from  12  A.M.  to  12  P.M. 


t 

24v 

h 

24 

X 

0.0469 

=  1.13 

24 

X 

0.2308 

r  5.54 

h 

24 

X 

0.5000 

=12.00 

24 

X 

0.7692 

=  18.46 

*5 

24 

X 

0.9530 

=22.87 

The  best  times  during  the  day  to  take  the  measurements  are  therefore 
1:08  A.M.,  5:32  A.M.,  12  H00N,  6:88  P.M.,  10:52  P.M. 


49. 

• 

CHAPTER  IV 
ACCURACY  OF  INTEGRATION  FORMULAS. 

19.  Introduction:     In  this  chapter  the  inherent  accuracy  of  the  integration 
formulas  is  to  be  discussed.    It  must  be  remembered,  however,  that  if  the 
given  ordinates  are  results  of  observations  they  are  liable  to  error  and  it  " 
is  of  little  value  to  have  excessive  accuracy  in  the  formula. 

20,  Relative  Accuracy  of  Trapezoidal  Rule,  Simpson's  Rule,  and  'Veddle's 
Rule.    In  order  to  compare  the  accuracy  of  the  three  formulas,  the  same 
interval  and  the  same  number  of  subintervals  must  be  taken  in  each  case. 
Six  is  the  smallest  number  of  subintervals  that  can  be  taken  so  the  interval 
is  chosen  as  6w.    m  is  taken  as  the  midpoint  and  the  points  of  the  sub- 
divisions will  be  x  z  m  -  3w,  m  -  2w,  m  -  w,  m,  m  4  w,  m  f  2w,  m  4  3w. 

Assume  that  the  function  y  s  f (x)  to  be    iifegrated  is  continuous  and 
has  continuous  derivatives  throughout  the  interval  and  that  the  value  of  w 
is  such  that  f (x)  can  be  represented  in  the  interval  by  a  convergent  power 
series  in  w  given  by  a  Taylor  expansion  of  the  function. 
fm  4  3w 

Let/          f(x)  dx  s  F(m  4  3w)  -  F(m  -  3w)  be  the  true  value  of  the 
Jm  -  3w 

integral. 

By  expanding  into  a  Taylor  series 

F(m  ±  3w)  s  F(m)  ±  3wf(m)  +  |  w2f'(m)  ±  .§7w3f"(m)  +  fiw4f-  (m) 

±  243wVv(m)  4  ££w6fv(m)  ±  ££Vf vi (m)  4- 
120                    720  5040 

Thus  the  true  value  of  the  integral  is 
Cm  4  3w 

(20:1)  /           fix)  dx  =  6w  f(m)  ♦  9w3  f»(m)  4  fiw5  fiv(m)  +  §|§w7  fvi(m)  f... 
J  m  -  3w 

By  applying  the  trapezoidal  rule 

fm  f  3f  (x)  ix  a  |  j~*(m  "  3w)  *  f  <m  4  3w)  4  2  £f  (m  -  2w)  4  f  (m  -  w)  *  f  (m) 
Jm  -  3w  2  ^* 

4  f  (m  4  w)  4  f  (m  4  2w))j 
If  f  (m  -  3w),  f  (m  4  5*),  f(m-2w)t...  are  replaced  by  Taylor  expansions,  then 


50, 


f m  4  3w  , 
(20:2)    /  f (x)  dx  5  6w  f (m)  4  -±|w°  f"(m)  + 

J  m  -  3w 


fiv(m)  4  ^£w7  fvi(m) 
8  120 


If  (20:2)  is  subtracted  from  (20:1)  then 


34„7  ^Vi 


•  •   •  • 


3y  applying  Simpson's  one-third  rule 

Cm  +  3w  r  r 

(20:3)  /  f (x)  di  :  J    f (m  •  3w)  ¥  f (m  +  3w)  4  4    f (m  -  2w)  4  f (m) 

J  m  -  3w  3  L  <• 

+  f(m  -  2ir)j    4  2  [f (m  -  w)  4  f(m  4  w)]J 
If  the  functions  on  the  right  are  replaced  by  Taylor  expansions  then 

(20:4)  /  m  +  3f(x)  dx  =  6w  f(m)  +  9W3  f"(m)  4  ||*5  fiv(m)  +  §§§w7  fvi(m)4.. 
J  m  -  3w  12 

If  (20:4)  is  subtracted  from  (20:1)  then 

S      z'Zt  fiY(m)  -  12*1  fVi(m) 
■l      30  420 

3 

By  applying  V'eddle's  rule 

/      +  3f  (x)  dx  =  |S  \f[m  -  3w)  4  f  (m  4  3w)  4  5  [f  (m  -  2w)  4  f  (m  4  2w)j 
Jm  -  3w 

4  f (m  -  w)  4  f (m  4  w)  4  6  f (m) 
If  the  functions  on  the  ritht  are  replaced  by  Taylor  expansions  then 

(20:6)     r  *  3?(I)  dx  a  6w  t  (m)  +  9w3  fw(m)  +  |lw5  fiv(m)  +  ?w?  fvi(m)+##- 
jm  -  3w  ^u  8 


51. 


If  (20:6)  is  subtracted  from  (20:1)  then 
(20:7)  =  -  -jJL  fvi(m) 


By  comparing  the  three  formulas 


20:3 
20:5 


„W?  f«(m)  -  11|  ^  fiT(B)  .  _34w7  fvi(m) 


8 1  " 

3 


40 


30 


420 


E  -Zllfvi(m) 
W      "  140 


20:7 

it  can  be  seer,  that  if  w  is  small  enough  so  that  the  Taylor  series  con- 
verges rapidly  ,       Weddle's  rule  is  the  most  accurate  and  Simpson's  rule 
more  accurate  than  the  trapezoidal  rule. 
21.    Formula  for  Error  in  Trapezoidal  Rule. 
Series  formula. 


Let   /  f  (x)  dx  -  F(m  4  w)  -  F(m)    be  the  true  value.    Expanding  the 
J  m 

right  side  by  Taylor's  theorem,  the  result  is 

(21:1)  I  :  I  f(x)Wdx  =  w  f(m)  +  |Jf ' (m)  +  ^M(m)  +  |Jf"  (a)  ♦  §^iv(«0 
J  m 

*|^(»)  *  ^Vi(m) 
By  applying  the  trapezoidal  rule 


(21:2)  y? (x)Wdx  s  |  [f(m)  t  f  (m  +  w)] 


Replacing  functions  on  the  right  by  Taylor  expansions  the  result  is 


(21:3)  Ir 


f  m  •»» 
=  /  f(x 
J  m 


x)  dx  =  w  f (m)  tlf'li)  +  2L.  f"(m)  +  JL  f"'  (m) 


12 


^flTtm^4fT(n)  +  4-ofT1(m)  +  --- 


(21:4)  I-J      Er^-^f^mjjslfm  (m)sJ^f i^(m)-lLfT(in) ^sl_fvl (m)- 

T       X»     12  24      V     "80       v   '  360     V   ;^20le       1  ' 

This  Is  the  error  for  one  aubinterval  from  x  =  m  to  x     m  »  w. 


To  get  the  error  for  the  whole  interval  from  x  r  x0  to  x  z  xQ  +  nw 
let  m  =  Xq,  x^,  X£,  •  •  •  •  *  Xn-l  8111(1  add  the  results. 


w3  r 


-  w4 
w£ 

8 


(21,5)  Thus  ET<>n  =  -  XT  l_ftt(xo)  *  f"(xl)  '"(^-i) 

>'(x0)  *  f  (xx)  *  .  .  .  ♦  fn'(*n-l) 
fiv(xQ)  a  fiv(Xl)  f^x^)] 


where  X  is  some  unknown  value  of  x  between  xQ  and 
Examples. 

1.    Compute  the  inherent  error  in  the  answer  to  (11, example  7) 
f(x)  =  ±,      f"(x)  =  4,      f-  (x)  -  =f 


1/513  -      .0000075  1/514  =  .0000001 

1/523  -      .0000071  1/524  ■  .0000001 

1/53 3  I      .0000067  1/53}  -  .0000001 

1/543  .      .0000063  1/544  I  .0000001 

ET?=  "  6  C0000276)  _  I  (,.0000004) 

4 

•  -.0000046  -  .0000001  =  -.0000045 


22«    Formulas  for  Error  in  Simpson's  one-third  rule. 


a*    Series  Formula  for  Error. 
Let  m  be  the  midpoint  of  interval  of  width  2w.    The  true  value  of  integral 


rni  m  fm-fw 

/         f(x)  dx     is  /         f(x)  dx  =  F(m  +  w)  -  F(m  -  w).  (22:1) 
J m  -  w  -  w 

Expanding  "by  Taylor's  theorem 

2!       —      3 J       *  *  4! 


F(m  ±w)  =  F(m)  ±  w  f  (m)  +  |_  f  •  (m)  ±21-  f»(m)  +       f*  (m)i  .  .  .  . 


Tm  "fw  ,  -j  5    <  7  « 

(22:2)  Thus  /         f(x)  dx  =  2  I  w  f(m)  +  TL  f(m)  f  |-  flv(m)  +       f V1  (m) 


Applying  Simpson's  one-third  rule 

(22:3)  J  m  +f *x)  dx  s  |  [f  (m  *  w)  *  4f  (m)  *  f  (m  ♦  w)  j 
Replacing  functions  on  right  by  Taylor  expansions  the  result  is 
(22:4)    (  m  *  f(x)  dx  -  ^  Uf  (m)  ♦  w2  f»(m)  +  |£  fiv(m)  t  |ji  fvi(m) 

If  (22:4)  is  subtracted  from  (22:2)  then 

(m.b)  v  ■  -  i  flv<B>  -  i§5o  fTi'»)  -  fVlii(m>  

3 

This  is  the  error  for  two  subintervals  from  x  :  1  •  f  to  1  *  n  f  v« 
Tbe  error  of  the  interval  from  xQ  to  xQ  +  nw  can  be  obtained  by  putting 
m  s  x^,  x  ,  x  ,  .  .  .  •  ,  x    ,  in  (22:5)  and  adding  the  results. 

(22:6)     Thus  E6i  =  -  ||  [flT(xl)  +  f1'!^)  *  •••  t  ^''Vl'] 

1890  [fTl,Xl»  +  fVl(X3'  *  "•  +fTl(lB-lJ 


"  90 
3  _7 


90720 

5 


(22:7)  or    EH    =  -  J2L  fiv(X)  -  Sill  fVi(X)  -  -iilL  fviii(X) 
81  9  1890  90720 

3 


where  X  is  some  unknown  value  between  xQ  and  xQ  +  nw. 
b.    Formula  In  Terms  of  Differences. 

If  the  analytical  form  of  the  function  is  unknown  or  if  it  is  difficult 
to  calculate  the  higher  derivatives  then  the  formulas  (22:6)  or  (22:7)  are 
useless  in  that  form.    However,  if  the  derivatives  are  replaced  by  differ- 
ences then  the  error  may  be  calculated  when  only  the  ordinates  are  given, 
remembering  of  course  that  there  are  probably  errors  in  the  ordinates  them- 
selves if  found  from  observation. 

To  replace  the  derivatives  by  differences  Stirling* s  interpolation 
formula  may  be  written  to  eighth  differences,  differentiated  with  respect 
to  x,  where  x  =  m  +•  uw  or  u  «  x  ^  mt  and  then  u  -  0  be  put  in  each 
derivative. 

Thus     fiT(M)   »  ^(Am.Sw  -  ^V3w  ♦  -g^Vto  '    •  •> 

w 

If  m  s  x^,  x3,  •  .  .  •  ,  xn-1  are  placed  in  formula  (22:6)  and  the 
eighth  differences  neglected,  then 


54. 


(22:8)  \  =  "  foU4y-l  +  *\  t  aS3  4 
3 

.6 


•  •  . 


+  d°y0  ♦  d6y2  *  .  .  .  +  A6yn_4) 
From  this  formula  for  error  it  cank  sWt  iiat  the  central  difference 
integration  formula  (12:1)  is  simply  Simpson's  rule  plus  correction  terms, 
c.    Formula  in  Terms  of  Given  Ordi nates. 

The  differences  in  (22:8)  may  be  replaced  by  their  values  in  terms  of 


y's.    It  is  of  little  value  to  keep  the  sixth  differences  because  the  y*s 
obtained  by  observation  are  liable  to  error  and  errors  become  greater  as 
the  differencing  is  continued. 

44y-l  "  ?3  ~  4y2  *  6yl  "  4yo  +  y-l 
A^l    =  yK  "  4yA  t  6y,  -  4yp  +  y 


55, 


y7  -  4y6  +  6y5  -  4y4  *  y3 


A  yn-3  ■  yn+i  -  4yn  *  6yn-i  *  *yn-2  +  yn-3 

Thus,  by  substituting  in  (22:8) 


I     =  "  |o  [y-l  +  yn*l  '  4(yo  +  yn}  +  7(yl  +  yn-l}  "  8(y2  +  y4  +  •  •  +  yn^ 
3 

*  8(y3  4  y5  t  .  .  ,  ♦  yn-3)j 
when  n  ^  6. 

E*l  5  "  lo  [y-l  4  y3  -  4(y2  +  yo}  4  fyl]      When  n  * 


E, 


lx  *  •  h  f  75  '  4(yo  +  y4}  +  ?(yl  +  y3}  *  8y2] 


when  n  a  4 


d»    Chevilliet's  Formula. 


If  in  formula  (22:2)  m  is  replaced  by  x^,  x»f  .  •  *  ,  xn^t  the 
intervals  taken  as  x0  to  Xg,  x2  to  x4,  •  •  •  ,  x^g  to  xn  and  the  results 
added  then 

(22:10)  I  s  J  nf(x)  dx  =  2w[f(x1)  +  f(x3)  *  4  *  ♦  ♦  ff*^)] 

*  4  lflT«xl>  +  '^V  4  '  '  •  +  fiV<*n-l>] 


The  first  line  represents  the  sum  of  the  areas  of  n/2  rectangles  with 
bases  of  width  2w  and  with  altitudes  ffx^),  f(x3),  •  .  •  ,  f(xn„i)»  Tne 


remaining  terras  are  corrective  terms. 


i: 


ft 


The  area  ABCDEFGHA  represents  the 
sum  of  the  areas  of  n/2  rectangles 
and  areas  HC I- IDN-N JE+  JFK-KGL-BHM 
are  represented  by  corrective  terms. 

The  formula  (22:10)  may  there- 
fore be  written  as  E  +  Cj^  +  Cg  +  ... 
where  R  stands  for  the  sum  of  areas 


of  rectangles  and  C^,  Cg,  •  •  •  are  the  corrective  terms. 
The  formula  (22:6)  may  be  written  as  follows 

E 


31  - 

3 


|[2w  f1^)  *  2w  fiv(x3)  *  .  .  .  t  2w  fiT(xn.1)] 
,  |*L  ^fviii(Xi)  +  fviii(x3)  +  .  .  .  +  f*"«W]*  .  .  .} 

The  terms  in  the  first  brar?  represent  sum  of  the  areas  of  n/2 
rectangles  with  altitudes  fiv(x,),  fiv(x3),  •  .  ,  fiv(xn«i)  and  with  bases 
of  width  2w.    The  terms  within    he  other  braces  are  correction  terms.  Thus 
all  the  terms  within  the  braces  represent  a  close  approximation  to  the  area 
under  curve  y  ■  fiv(x)  from  x    to  xn. 


( 22*11] 


The  error  is    [jR  t  C±  +  C2  t  .  .J  - 


7  1      42  d 


...  • 


or 


7  x  42 

If  w  is  small  enough  so  that       and  Cg  are  small  then  the  area  under 
the  graph  will  be  closely  approximated  by  £f 


7^1  T  4T2 

The  true  value  of  the  area  under  curve  of  y  =  f    (x)  is 


/   Vv(x)  dx  =  [f'"(x)Tn  s  f  (x  )  -  f™(x  )  • 

Substituting  this  value  in  (22:11)  the  result  is 

(22:12)  -*l_[f-(xn)  -  f-(xQ)J 

3 

This  is  known  as  Chevilliet*s  formula.    It  is  unreliable  unless  w  is  small, 
but  it  is  useful  in  finding  the  value  of  w  for  any  desired  degree  of 
accuracy. 

22.    Example  1. 

Compute  error  in  answer  to  (11  example  2)  by  Series  formula  (22:6) 

fiT(x)  =  m 

X6 


2 


 1    .030  +  .003     =  .000  +  .    Thus  there  is  no  error  in  (11  example 

*1       90  L  J 


2b)  to  three  decimal  places. 
Example  2. 

Compute  error  in  answer  to  (11  example  4)  by  (22  :B)* 


E_    =  -  _J_(19  -  14  -  8)  +  -i — 

540  4536 


3 


(-  59  +  46  4  27) 


57. 


=  <+  .007  -  .003  =  .004      There  is  no  error  in  the  first  decimal  place. 
Example  3.      Compute  error  in  (11  example  2)  by  (22:12) 

E  --— ' -24-24  )  r  24  (.00002-  .00411)  -  2  (-.00409)  «  -.00054 

BA  -  iec  -98  w    im  T5" 

3 

There  is  no  error  in  the  answer  to  three  decimal  places. 

Example  4.    Compute  error  in  the  answer  to  (11  example  7)  by  (22:9) 


*  Refer  to  page  30,  example  2,  for  the  table  of  differences. 


.020000000  s 
.017857143  = 
7(.019230769)s 
7 (.018518519 )« 


. 020000000 
.017857143 
.134615383 
.129629633 


f  .302102159 


-4 (.019607843)= 
- 4 (. 018181818 )s 
-8 (.018867925)= 


-  .078431372 

-  .072727272 

-  . 150943400 

-  .302102044 


.302102159 
-  .502102044 

A-  (.000000115)  =  .000000001277 
90   

e.    Formulas  in  Terms  of  Twc  Computed  Results. 

If  two  computations  of  a  definite  integral  are  made  by  Simpson's  rule 
but  using  different  values  for  w  a  comparison  of  the  errors  in  the  two 
results  can  be  made  by  Chevilliet's  formula. 

Let  R^  and  Rg  be  the  two  results  obtained  when       and  wg  are  used.  The 
errors  are       and  Eg  respectively. 


Then  by  (21:18) 


(22:14)  Thus  if  w?  s         then  E,  =  iA 


n 


Let  I  be  the  true  value  of  the  given  integrand  and  then 


(1)     I  =  R1  +  Ej_  =  Rx  +  m^g 


(2)     I  s  Rg  +  Eg 


Subtracting  (1)  from  (2) 


Thus,  if  we  take  4  times  as  many  subdivisions,  that  is  m  =  4,  the 


3  1  i 

result  will  be    Eg  s    or  E2  will  be  about  -z=-th  of  the  difference 

255  255 


of  the  two  results. 
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These  formulas  are  accurate  only  if  Chevilliet's  formula  is  accurate 
and  this  depends  upon  the  value  of  w. 

f.    To  Find  Value  of  w  for  a  Specified  Degree  of  Accuracy  in  the 
Integral.  This  may  be  found  by  ChevillieVs  formula.    If  the  analytic  form 
of  the  function  is  known  and  the  third  derivative  can  be  found  then  w  may 
be  found  by  substituting  in  formula  (22:12)  and  solving  for  w.    If  the  ana- 
lytical form  of  function  is  unknown  or  if  the  third  derivative  cannot  easily 
be  determined  then  it  is  necessary  to  assume  a  value  w^  for  w  and  to  find 
E^  by  formula  (22:9).    Then  if  2^  is  the  specified  error  in  the  computed 
integral  by  (22:13) 

Wl  ^ 

(22:16)        jj-  =  -5-    and  w  :  w1  —J 

w  and       must  be  small  enough  to  give  an  accurate  result  in  Chevilliet's 
formula. 

(22)    Example  1.    ?/hat  is  the  velue  of  w  necessary  to  give  the  integral 

I     —    (11  example  2)  correct  to  5  decimal  places, 
73  x2 

The  error  b\   (,2  2:9)   is  ♦00040644  = 

E„  K  .000005  4. 

43                                         m  ,  w. 000005,* 
9       1                                  <  1(.000406}    8  *Si 

As  x0  +  nw  -  x0  *  nw      we  must  divide  the  interval  (3,9)  into  more 

than  18  subdivisions  to  get  the  result  correct  to  five  decimal  places. 

23.    Formula  for  Error  in  'Veddle's  Rule.    By  formula  {10:7) 
7 

E_  »  -     —  fyi(m)    for  one  subinterval  from  x  =  m-3wtox»m  +  3w. 
w  140 

To  get  the  error  of  the  whole  interval  from  xQ  to  x_  put  m  =  x, ,  x  ,  ,  ,  , 
xn-3* 
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Then  Sw  -  -         ffvi(x5)  ♦  fti^)  ♦  .  .  .  ♦  fVi(x^)] 


(25:1) 


or  - 


'rr7 


7  „vi 


fvl(X) 


where  X  i9  some  unknown  value  between  x0  and  xn. 

Example.      Compute  error  in  answer  to  (11  example  2)  by  (23:1) 

fvi(x)  m  ^040 

x 

24.    Remainder  Terms  in  Central-Difference  Formulas.     The  remainder  term 
in  Bessel's  interpolation  formula  is 

w2n  +  2f2n  +  2^      -  _         2      g  InT-l2 
Rn  (2n*2)!     *     ^         *A         V  *  '  *  CV   5  > 

As  the  central-difference  integration  formulas  was  obtained  by  integrating 
f (x)  dx  the  error  in  the  formula  may  be  obtained  by  integrating  Rn(x)  dx. 
As  x  «  x0+"w(v*i)  in  Bessel's  interpolation  formula  then  dx  -  w  dv  and  the 
jrror  for  a  single  subinterval  of  width  w  is 


(24:1) 


)    A0  +  w 

E  -  /     R^x)  dx 
J  x0 


*  ..2n  +  2*2n  +  2 


(2n  +  2)1 


w—  -  -f 
(2n 

(v2  _  sn 


^  (v2  -i)(v2  -|)  .  .  . 


2 


2w2n  ;  3f2n  ♦  g^gj/       2  _         2  _  p 

(2n  +  2}  I V  4A         J)  •  •  • 


(24:2)     If  V 


.A 


(T2  .  2n_pf}  dv 
/    (v2  -  i)(v2  -  |)  .  .  .  (va  . 

'f2n;     |v  i 


— 5 — )  dv 


w2n  +  5f 2n  +  2 
(2n  +  2) 


Then  E 

This  is  the  error  for  single  interval  of  width  w. 

Let  M    denote  the  maximum  value  of  f  (x)  in  the  interval  (x0,  xn). 


Then  as  there  are  xn  ~  xo    subintervals  from  x  -  x0  to  x  «  xn  the  total 

w 

error  is 

If  all  differences  higher  than  the  fourth  in  formula  (12:2)  are 
neglected  and  n  made  equal  to  2  then  from  (24:2) 


'    »  /  (v2.-  i)(v2  -  £)(v2  .  25)  dv  .  _  1£1 

"Jo  *  V  16S 


Thus  by  (24:^) 


n 


(24:4)     ^  EXo<12(x2-x0) 


(24:5)    or  e£<  (xa  -  x0) 

where  4^y  is  the  largest  of  the  sixth  differences  in  interval  (x0,xn) 

If  all  differences  higher  than  the  sixth  in  (12:2)  are  neglected  then 
n  ■  5  and  (24:2)  becomes 


Substituting  in  (24:5) 


Exn 


2497wSM: 


5626600 


(Xo    -  XQ) 


where  A°y  is  the  largest  of  the  eighth  differences  in  the  interval  (x0,xn). 

When  all  differences  higher  than  the  fourth  in  formula  (12:1)  or  the 
third  in  formula  (12:5)  are  neglected  then  formula  (24:5)  is  to  be  used  in 
finding  the  error.    When  differences  higher  than  sixth  are  neglected  in 


(12:1)  or  the  fifth  in  formula  (12:5)  then  formula  (24:7)  is  to  be  used  in 
determining  the  error. 

These  are  the  formulas  to  be  used  in  computing  error  in  the  integration 
formula  (12:5)  derived  from  Bessel's  interpolation  formula.    They  may  also 
be  used  in  computing  the  error  in  the  integration  formula  (12:1)  derived 
from  Stirling's  interpolation  formula  since  (12:1)  is  at  least  as  accurate 


24.    Example  1.    Compute  the  inherent  error  in  the  answer  to  (12  example  2). 


as  (12:5). 


25.    Remainder  in  the  Euler-Maclaurin  Formula. 


then  by  putting  z 


■  a  +  w 


u 


By  expanding  into  a  Taylor  series  and  expressing  the  remainder  as  a 
definite  integral 


By  replacing,  in  the  above  formula,  #?by  <p',  .  .  .    (p  2p_1 


and  at  the  same  time  replacing  2p  by  2p  -  1,  2p  -  2, 


then  it  follows 


Jordan1 8 

Cours  D'Analyse",  vol.  II  pg.  100  Published  I885. 
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1 

that 

<p  (a  +  w)  -/(a)  =  w  <^(a)  +  ^  ^*(a)  +  

<—  • 

/(a  ♦  w)  -9?(a)  -  w^(a)  ♦  |£  <^a)  ♦  ....+  »2g~2  ^2p(a) 

/0>1(a  •  ■  -  u>  iSri-  * 

<^2P"1(a  +  w)  -<^2p*1(a)  =  w^2p(a)  +  J  ^2P+1(a  +  w  -  u)  u  du 

By  adding  these  equations  to  equation  (2^:1)  after  having  multiplied 

them  respectively  by  their  coefficients  A^w,  A£W2,   ,  ■^2p-lv'^^"^, 

chosen  in  such  a  way  that  the  terms  <^(a) ,  Cp  '(a) , . . . .  ^>2p(a)    disappear  in 
the  right  member  of  the  equation. 

(2^:2)^ (a  +  w)  -^(a)  +  A1w^'(a  +  w)  -  (^'(a)]  +  A2w2|V*(a  +  w)  -<?*(a)]  + 
 ♦  A^p-!  w2P"3[^2P-1(a  ♦  w)  -^^(a)]    =  w/(a) 

+^^2p+^(a  +  w  -  u)  F(u)  du,  where 

f(u)-ll?TTtAi>'715iITT^-"+A2p.1w2P-1u  ■ 

The  values  of  the  A's  are  determined  by  the  following  equations:* 

1.  +  Ai   -  0 

2.  1 

1  Al 

1  Al 
(2p7T  +  T2FTTT +  •  •  •  •  +  A2P-i  »  0 

*  "Cours  D'Analyse",  vol.  II,  by  M.  C.  Jordan,  pg.  101  -  Published  188?. 

which  give 

A1  ■  -L  A 2-  -  0, 


(-l)^"~^B 

•  A2p-1  "  0.  A2  =  =    ~2k),  k 


where  the  B's  are  the  Bernoullian  numbers. 

Let  f(x)  be  equal  to  the  derivative  of  ^(x),  Then 
ra+w 

(a  +  w)  m  I     f(x)  dx.    Substituting  this  and  the  values  of  the 

A's  in  (2^:2)  the  following  is  given: 

(250)y  &f(x)  dx  =  w  f (a  *  w)  ♦  f(a)  _  |  B^  fI(a  +  w)  _  f ,(a)j  +   

+  (-l)P'l3p-lV/2?"2[f2p^(a  +  w)  .  f2p-5(a)]  +  d 
(2p  -  2)1         I  J 


where 


It  can  be  seen  from  the  following  development  that 


(25:5) 


!    -  "r-i  1  w2u2F-2  ? 

T^FTT  ~  ff  ;2P-d:     2:  (2F-2)i   J 


is  the  same  as 


(25=5)      W"  <p^). 

<^(n),  £^(n),  ...  are  defined  as  the  coefficients  of  the  powers  of  x  in 
.    nx  _  ■> 

the  expansion  of    and  are  called  3ernoullian  polynomials. 

ex  -  1  n2x2  J.J 

nx  +  —r-  +  — -t  ... 


(25:7) 


-  1      enx  -  1 


(1  -  *+  *£. 


-  1 


4 


+  •  •  •  • ) 

■  n+  </|(n)X>  <^n)x2  +  (^(n^x"  + 

ottl      1     •*  Bin04"' 

where*/  (n)  -  ,    •  »    ,  -  ±  &■  4  ^  r  -   


If*-  2p  -  1  and  n  -  -£    then"/1,    A-)  -      u2p  1  «-!  

Thus  (2^:5)  is  the  same  as  (25:6)  and  (25:4)  can  be  written  as 

(25:3)      R  - y*f2p<a  +  w  -  u)  du  w^^.^H)  du. 

In  (25:6)  u  varies  from  0  to  w  and  H  varies  from  0  to  1, 

w 

If  ju  is  taken  as  an  intermediate  value  between  the  maximum  and  minimum  valuee 
of  the  factor  f 2P(a  +  w  -  u)  then,  as  the  argument  (a  +  w  -  u)  varies  from 
(a  +  w)  to  (a)  in  the  course  of  integration,  y.  *  f^p(a  +  -  i)  where  Q  ia  a 
value  between  0  and  1  and  R  can  be  written  as  follows: 

R  -   w2V2p-i^)  du- 

Then  if  t  is  substituted  for 


3y  taking  the  derivative  of  (25:7)  in  respect  to  n 


xe 


nx 


e    -  1 


1  +  <p[(n)x  +<f'0(n)  x2  ♦  .  .  .  +^'(n) 


Al30, 


xe1^    ^  - 

iX  -    1  *     L:x  -  1 


-   1  X 


ex  -  1 


-  x  [n  ♦  <^(n)x  +  .  .  .  +<£(n)x°+  ■  •  •  ]  +  [  ~  f  4  1 


2  A 
BiX  3oX 


i-iX"      d-x  7 

-  t —  +  .  .  .  / 

21      hi  J 

Upon  equating  the  coefficients  of  the  same  powers  of  x, 


(2k)i 


Thus    wl    cp^it)  dt  «    J  jf&M  + 


(-DP  Bs 
(2p)J 


J 


dt 


It  can  be  seen  from  (2^:7)  that  ^(n)  vanishes  for  n  »=  0  and  for  n»l, 
enx  _  i 

for  in  each  case  ±    =  n. 

ex  -  1 

Then^2p("0  vanishes  at  the  limits  0  and  1  and 

R-f*(aW-1)P  BP  =5^ 

(2p)J 

To  apply  this  formula  to  /    f(x)  dx  it  is  only  necessary  to  place 

J  a 

N  *  b  -  a.    To  get  more  accurate  results  the  interval  from  a  to  b  may  be 


divided  into  subintervals.    Then  w 


b  -  a 

n 


and 


/  b  /"a+w  ra+2w  ra+nw 

/      f(x)  dx  -/     f(x)  dx  +   /     f(x)  dx  +  .  .  .  +  /  f(x)  dx 

'a  7a  7a+w  ,/a+(n-l)w 

If  formula  (25:3)  is  applied  to  each  of  these  integrals  and  the  results 

added  then 

]tM  dx  .  „'[*  f(a)  ♦  f(a  4  w)  *  f(a  .  — 1  w)  +  i  f(b)j 

,     ,  r  7  H**"1  3p-i  w2p"2 

■iB^fW-fW/  — 

[f2P-5(b)-f2I"-2(a)] 
This  is  the  Euler-Maclaurin  formula  and  the  remainder  is  given  by 
(_1)P  3n  w2P+1 

(25:9)    R  =  au   i)1   where  p.  is  some  value  between  the 

maximum  and  minimum  value  of  f^P(x)  in  the  interval  from  a  to  b« 

The  simplest  and  most  useful  rule  for  obtaining  the  most  accurate  result 
from  the  Euler-Maclaurin  formula  has  been  given  by  Charlier*  as  follows: 


♦"numerical  Mathematical  Analysis"  by  Scarborough,  pg.  169. 


"In  stopping  with  any  term  in  Euler's  formula  the  error  committed  is  less 
than  twice  the  first  neglected  term"    Thus  the  most  accurate  result  is 
obtained  by  stopping  with  the  term  just  before  the  smallest  term. - 
26.    Formulas  for  Error  in  Gauss'  Integration  Formula.  As  the  derivation  of 
the  formulas  are  long  they  have  been  omitted  and  only  the  formulas  them- 
selves given. 

(26:1)    a.*  En  -   2—  [   al  .  J2 

*  Q     2a  +  1  I  l-?-=  .  .  .(2n  -  1)J 


where  the  b's  are  the  coefficients  in  the  power  series  f(x)  ■  b0  +  b-^x 
o 

+  +  •  .  .    and  the  values  of  pc,  jS^>,  3^  .  .  .  are  as  follows: 

/q       n(n  -  1)  a       n(n  -  1' 
/®2  "  2(2n  -lj  Po  '  2(2^ 


fl       n(n  -  1)     _       n(n  -  l)(n  -  2)(n  -  3)  * 

«  "2(n  ~  l)2       n(n  -  l)(n  -  2)(n  -  3) 
22(2n  -  l)2      2-4(2n  -  l)(2n  -  5) 

This  is  used  when  the  limits  of  integration  in  Gauss'  formula  have  been 
taken  as  -1  and  1. 


*  Derived  in  Todhunter's  "Functions  of  Laplace,  Lame,  and  Bessel"  pg.  106. 


(26:2)      b.*    EL  =   -b 

(2n 


♦  l)22n  (  1-5-5  ...(2n  -  1)]  + 


-    :  : 


L 


(n  +  l)(n  +  2)      n(n  -  1) 
+    2n  -  1 


2n  +  5 


where  the  b's  are  coefficients  in  the  power  aeries 
Cf(u  +  |)  =  bQ  +  b,w  +  b2u2  +  .  .  .  +  ^2nu2n  +  .  .  .  . 

This  is  used  when  the  limits  of  integration  in  Gauss  1  formula  have  been 
taken  as  0  and  1. 


(26:5)      c"Er  = 


r*  v2n+l 
(b  ~  a) 


G        (2n)l  [in  +  l)(n  +  2;   ...  2nJ     2n  +  1 
where  £  is  some  value  of  x  between  a  and  b.    f (x)  must  not  change  sign  in 
the  interval  of  integration  (a,b). 

(26:4)    d***    S    -  f/Kffi  *       2     [,      5j  J2 
^        ^  G        (2n)  i      2n  +  1  |_(2n  -1)1. 

where  £  lies  between  -1  and  1. 

26.  Examples: 

1.    Compute  the  inherent  error  in  the  answer  to  (18  example  1)  using 
(26:1) 

f(x)  -  §<1  ♦  f)"1  .  |[l  -  (§)*(*)2-(§>*  ♦  (f)4  -  (|)5  ♦  (f)6  - 

♦  Cf)8j 
b4  "  I  <:i->-  I  (15) 


*  Derived  in  Todhunter's  "Functions  of  Laplace,  Lame,  and  Bessel.  pg.  103. 
**  Given  in  Scarborough's  "Numerical  Mathematical  Analysis."  pg.  I67. 
***  Goursat-Hedrick  "Mathematical  Analysis"  Vol.  I.  pg.  207,  ex.  24. 


B„-l.      B2-i,  8,-1 
When  substituted  in  (26:1) 

E    =  W  (1)2/  1  +    Xfl^  f  ll    1   f?'^-6  ♦  134  ♦  I 

=  (.088888)  I  .0625000  t  (.015625) (1.190475)  +  (.005906) (1.111040) J 

-  (.088888) (.08 1505)  =  .0072446 
Example  2.    Compute  the  inherent  error  in  the  answer  to  (18  example  1)  using 
(26:4) 

n  =  2  fiv(x)  -   ^  

(x  ♦  2)5 

fiv(f)j  where  €  -->!-  24 

fiv(£).  where  J  ml.  mM  m^m  .O9876 

when  I  «  -1S  f  [^]2  -  -044444 


1  2(2 


.0001 


10955 

The  error  lies  between  ,0444  and  .0001.    It  can  be  seen  that  the  error 
cannot  be  found  very  accurately  by  this  method. 
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27.    Conclusions.    If  the  analytic  form  of  the  function  is  known  and  if  it 
is  easily  differentiated  it  is  advisable  to  use  the  Euler-Maclaurin  formula 
for  it  is  more  accurate  than  most  of  the  others.    If  the  analytic  form  of 
the  function  is  unknown  or  not  easily  differentiated  then  Gregory's  formula 
may  be  used. 

It  is  not  necessary  to  use  as  many  ordinates  for  the  same  degree  of 
accuracy  in  Gauss1  formula  as  in  the  other  formulas  but  the  labor  is  greater 
and  a  calculating  machine  should  be  used  for  exact  results. 

Newton-Cotes'  formula  and  the  trapezoidal  rule,  Simpson's  one-third 
rule  and  V/eddle's  rule,  which  are  special  cases  of  Newton-Cotes'  formula, 
are  less  accurate  but  require  much  less  labor.    They  are  accurate  to  a  few 
decimal  places.    Weddle's  rule  is  the  most  accurate  of  the  three  but  the 
number  of  ordinates  must  be  a  multiple  of  six  which  is  often  inconvenient. 
Simpson's  rule  is  nearly  as  accurate  and  can  be  used  for  any  even  number  of 
ordinates.    The  trapezoidal  formula  is  the  simplest  to  use  and  is  accurate 
to  about  two  decimal  places. 

The  central-difference  formula  derived  from  Stirling's  formula  (12:1) 
is  slightly  more  accurate  than  that  derived  from  Bessel's  formula  (12:5)  for 
it  converges  more  rapidly  and  it  uses  fewer  ordinate3  outside  the  range  of 
integration.    A  disadvantage  is  that  the  subintervals  must  be  even  in  (12:1) 
They  may  be  even  or  odd  in  (12:5).    Formula  (12:1)  is  Simpson's  one-third 
rule  with  correction  terms  and  (12:5)  i8  the  trapezoidal  rule  with 
correction  terms. 

The  accuracy    of  the  formulas  depends  much  upon  the  value  of  w  and  upon 
the  number  of  ordinates  used.    It  also  depends  upon  the  accuracy  of  any 
ordinates  derived  by  observation. 
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Appendix  A 
1«    Interpolation  Formulas. 

(7:1)    y  =  f(u)  -  y0  *  xy0  *            X)  &\           "  *>(f  ~  2>  &h0  +  •  •  • 

Newton's  Interpolation  Formula. 

(8:2)    y  =  f(a  ♦  xw)  =  yQ  ♦  x(A*-l  +^0)    +  ^ 

x(x2-l)      A-2                .    x2(x2  "  l)A-2 
51                  2  4J 

t  x(x2-  D(x2-                *A-2  t 

,  x(x2  -  i)(x2  -  4)(x2  -  9)  .  .      _  rn2  /^^/^^(n-i; 

• 

(2n  -  1):  2 
,  x2(x2  -  l)(x2  -  4)(x2  -  9)  ...  (x2  -  n  -  I2)  or, 

Stirling's  Formula. 

(3:4)    y  .  yG  *  y,  ♦  (x  -  *».  *  x(>'^  1}  &1  +o  ^2 

,  x(x  -  l)(x  -  t)  A5          x(x  -  lKx  ♦  l)(x  -  2)  A-2  *-l 

,  x(x  -  l)(x  ♦  lHx  -  2)(x  -  ft  £v  , 
5  J  —2 

x(x  -  lKx  +  l)(x  -  2)(x  *  2)....(x  -  n)(x+n-])  Z^V-n  + 

2n '  2 

+  (x  -  fe)(xKx  -  l)(x  -  2Kx  +  2)..(x  -  n)(x  +  n  -  l)  2n+l„ 

(2n+l)J                             *  Y'n 

Bessel ' s. Formula. 

(8:5)    y  -  70  !  71  +  v4yn  * 


(v2-i)  /y.i  */y0    v(v2  -i)  3 

a  J-i 


21 


+  4] — a  g    5]  ^^y-2 +  •••• 

where  x  ■  v  +  ■§-  Bessel's  Formula. 

Remainder  in  Bessel's  Formula 


n 


(9:3)    f(x)  =  ^  ■   'M^tf fftj 


where  f(xk)  -  Ak(xk  -  xQ)(xk  -  Xj)  ...  (xk  -  xk-1)(xk  -  xk+1)  ...  (xk  -  xn) 
and  F(x)  -.(x  -  xQ)(x  -  Jtj)  ...  (x  -  xn) 

Lagrange  1 s  Formula. 


72, 


2.      2.    Integration  Formulas. 


yf  x  +nw  r  o  ^      9      2.  4  3 

°y  dx  -  w  ny0  2^y0  ♦  ( -A  f )  AJfft  +  .  n5  +  n2)  4Z° 
xQ  t-        4  7  ■ 


^       2      +    5         5n  ;    hi  J 


Newton's  Formula. 


(11 


A0+nw 
:2)  /  yd: 


y  dx  =  |(yQ  +  2yx  ♦  2y2  +  .  .  .  +  2yn-1  +  f  ) 


Trapezoidal  Rule. 


f x0+nw 

(11:5)   /      y  dx  =  -  I  yQ  +  4(yx  +  y^  +  .  .  t  yn_i)  +  2(y2  +  y4  +  .  .  +  yn_2) 
J  xQ 


Simpson's  One-Third  Rule 


/x0+nw  r 
y  dx  =  K)  [y°  *  5yl  +  y2  +  6y5  +  y4  +  5y5  +  2y6  +  5y7  +  y8 

 + 

5yn-5  +  yn-4  +  6yn-5  +  yn-2  +  +  yn] 


+  6y9  +  y10  +  5yll  +  2y12  + 


x0+nw 


Veddle's  Rule. 


(12:1)  ^    >  dx  •  2w|yx  ♦  y5  ♦  y5  +  •  .  .  +  yn-1  +  ^(^y0  +4%  +  4% 

+  •  •  •  +  42yn_2)  -  lio^Vi  +  +  A5  +  •  +An.5) 
+  I5l2(A6y-2  +  *6yo  +  ^6y 


+  .  .  .  +  d6yn_4) j 


where  n  is  an  even  number.    Central  Difference  Formula  derived  from 

Stirling's  Formula. 


/  x0+nw 
(12:5)     /     y  dx  =  w 


n 


|2  ♦  yx  ♦  y2  ♦  .  .  .  +  yn.x  +  -f 


"  12*  2  '  2  ' 

U  /yn-2  +  dX-i  *  A^'-l 
♦   *) 

720  2  2 


191  rj^a^  !  ^"yn-2  ?  4^-2,  ] 

~  7M2?^  2  "  2  J  J 


where  n  is  an  even  or  odd  number.    Central  Difference  Formula  derived  from 

Bessel's  Formula. 


(13*4) 


y  dx 


J|2  +  yi  +  y2  +  .  . 
+  T20(yn-yo>  -  Toio(yn-yo)  +  -  -  •  j 


Euler-Maclaurin  Formula. 


(14:2)   /    y  dx 


=  N 


§S  *  ri  ♦  y£  ♦  .  .  .  *  yn_1  ♦  £2  -  yl^y^!  -  4y0) 


Gregory's  Formula. 
(15:2)  yQ  ♦  y^  *  7^  *  i      •    *  yn  -  m(yQ  ♦  yx  +  .  .  .  *  yj 


B  -  1 


2 

(y     +  y  )  -  m    ~  1  .  _  i\ 


-4  _  ,  .3 


720    m5Kyn  yoJ 


V/oolhouae's  Summation  Formula. 


(16=1)  y0  *  y1/m  ♦  y2/m  *  .  .  .  ♦  yn  .  m(y0  ♦  ^  ♦  ,  j.  ,  ♦ 

-  ±fMr.  ♦  yn)  -  -  -  4sT(^n-2  «  *%> 

-  (»2  -  i)(i>2  -  i)u?       .  ^  ) 

720  m'  ' 


480  m> 


Lubbock's  Formula. 


;;or"" 


(17:2)   /     f(x)  dx  =  K0f(x0)  +  H1f(xQ  +  w)  +  H2f(xQ  +  2w)  +  .  .  .  +Knf(x0+nv 


where 


EL*  -  (-l)n"k  w/  ^(2  -  1)  ...(2  -  k  +  1) 
'  Kl  (n-k)!  ^0 

(z-k-1)  ...  (z-n)dz 

Newton-Cotes  Formula. 


(18:1)  J    f(x)  dx  -  H0f(h0)  +  H^Ch^  .  .  .  Hf(ha) 


where  H,  * 


F(x) 

k     F'hJ^  x-hk 


dx 


and  h^  is  a  root  of  the  Legendre  polynomial  of  degree  n  +  1. 

Gauss'  Integration  Formula, 

Formulas  for  Errors  in  Integration  Formulas. 

(21:5)    Eyn  ~  "  f|  [*"(*©)  +  *"<*l)  *  ...  -  t9(ij} 

-  f£  [f  »'  (x0)  4  f  <»'  (Xl)  +  .  .  .  +  f«  (xjj 

-  f|  [f iv(x0)  -  f iv(xx)  +  *  f iv(xj] 


Error  in  Trapezoidal  Formula. 


(22:6)  Ef 


1/5 


-  ^[fiV(Xl)  ♦  fiV(x5)   ♦  .        .  4  ^(X^i)/ 

1890 


90720 


[fVi(x,)  4  fVi(x.)  ♦  .   .  .  +  fVi(xn-1)] 
[fVlii(Xl)   ♦  fViil(x5)   ♦  ...  4  fViii(xn^)] 


Series  Formula  for  error  in  Simpson's  One-Third  Rule. 
"1/5       90  1  —  -  r5  *•  *  '  '  +^Vj)J 


(22=8;    E„.  ._  -  -  |L  py.!  +  Ax  «  ^ 


Difference  Formula  for  error  in  Simpson's  One-Third  Rule. 


(22: 9)  \/*>  Hfe  [J-l  +  y*+l "  A(y° +  y»}  +  7(y! +  y"-1} 


'1/5 

-s(y2  +  + 

when  n  ■  6 


•  +  yn-2)  +  8(y5  +  y5  +  •  •  •  +  yn.O j 


when  n  »  2 


"  |o  [y-i +  y5  *  4(yo +  y4>  +  7(7! +  y5)  -  ey2J 


"V5 

when  n  ■  A. 


Formula  in  terms  of  given  ordinates  for  error  in  Simpson's  One-Third  Rule 
(22:12)     ^       -  -  £g[f  («,}  -  * *C*b>jr 


Chevilliet 1 s  Formula  for  error  in  Simpson's  One-Third  Rule. 
*2  "  Rl 

Formula  in  terms  of  two  computed  results. 


(22:15)    E2  » 


4 


m    -  i 


(22:16)     w  -  w. 


■4 


To  find  value  of  w  for  specified  degree  of 
accuracy. 


(25:1)    ^  .  -  3*    [  f»i(x?)  ♦  fi(x4)  ♦  ,  .  .  ♦  t'H^)] 


(24.5)  E^<W( 


x„  - 


516  %  n 


18. 


Formula  for  error  in  V/eddle's  Rule. 

xQ)    where  4^y  il  largest  of  the  sixth 
differences  in  interval  (x0,xn) 


Exo  ^    l4^Xn  -  x0)    where  ASyis  the  largest  of  the  eighth 

differences  in  interval  (x0,xn) 
Remainder  Terms  in  Central  Difference  Formulas. 


(26:1)    Er  = 


where  b's  are  coefficients  in  power  series  f(x)  «  bQ  +  b^x  +  b2x2  +  •  . 
and  fi's  are  as  follows  :  /30-/ 


n(n  -  1 
2      2(2n  -  1 


2(2n  -  if2      2  A(2n  -  l)(2n  -  5)  ro 


Formula  for  Error  in  Gauss  1  Formula. 


(26:2)  Er 


b  -  a 


ni 


(2n  +  l)22n(  l-5^2n  -  1) 


>2n*^ 


+  n(n  -  1) 

-    -  J 


(n  +  l)(n  *  2) 
2n  +  5 


where  b's  are  coefficients  in  power  series 

^(u  +  2)  ■  b0  +  b^u  +  b2u2  +  .  .  .  . 


Formula  for  error  in  Gauss  '  Formula 

2n+l 


(26.5)    E a-«£fDf,  ^—  I2  (*-») 

G      (2n)i   (,  (n  +  l)(n  +  2)  ,.2nJ      2n  +  1 

where  £  is  some  value  between  a  and  b. 

Formula  for  error  in  Gauss'  Formula. 


where      lies  between  -  1  and  i. 

Formula  for  error  in  Gauss'  Formula. 


Appendix  B.    Table  for  Values  in  Gauss'  Integration  Formula.* 


n  -  2 

vQ  =  C. 11270166 
v2  »  0.88729355 

1  9 

n  »  y 

VQ   »  O.O6945I84 

vx  «  0.55000948 
v2  =  O.66999052 
=  0.95056816 

R0  «       =  0.17592742 
R1  =.  R2  -  C. 52607258 

n  -  4 

v0  »  0.046910CS 
»  0.25076554 

v2  -  0.5 

=  0.76925466 

»  0.95505992 

R0  -  R4  -  0.11846544 
Rx  =  R^  =  0.25951454 
R2  *  0.28444444 

n  =  5 

v0  »  0.05576524 

V]>    «  0.16959551 

v2  =  C. 58069041 
=  0.61950959 
=  0.85060469 
=  0.96625476 

• 

R0  =  R^  =  0.08566225 
Rl  ■  R4  =  0.13058079 
R2  -  R^  =  0.25595697 

*  "The  Mathematics  of  Physics  and  Chemistry"  by  H.  Margeneau  and 
G.  Murphy.  Table  6,  pg.  464. 
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